CONTROL OF 
DISTRIBUTED 
PARAMETER AND 
STOCHASTIC 
SYSTEMS 



Edited by 

Shuping Chen 
Xunjing Li 
Jiongmin Yong 
Xun Yu Zhou 






SPRINGER SCIENCE+ 
BUSINESS MEDIA, LLC 




Control of Distributed Parameter 
and Stochastic Systems 




IFIP - The International Federation for Information Processing 



IFIP was founded in 1960 under the auspices of UNESCO, following the First World 
Computer Congress held in Paris the previous year. An umbrella organization for societies 
working in information processing, IFIP's aim is two-fold: to support information 
processing within its member countries and to encourage technology transfer to developing 
nations. As its mission statement clearly states, 

IFIP's mission is to be the leading, truly international, apolitical organization which 
encourages and assists in the development, exploitation and application of information 
technology for the benefit of all people. 

IFIP is a non-profitmaking organization, run almost solely by 2500 volunteers. It operates 
through a number of technical committees, which organize events and publications. IFIP's 
events range from an international congress to local seminars, but the most important are: 

• The IFIP World Computer Congress, held every second year; 

• open conferences; 

• working conferences. 

The flagship event is the IFIP World Computer Congress, at which both invited and 
contributed papers are presented. Contributed papers are rigorously refereed and the 
rejection rate is high. 

As with the Congress, participation in the open conferences is open to all and papers may 
be invited or submitted. Again, submitted papers are stringently refereed. 

The working conferences are structured differently. They are usually run by a working group 
and attendance is small and by invitation only. Their purpose is to create an atmosphere 
conducive to innovation and development. Refereeing is less rigorous and papers are 
subjected to extensive group discussion. 

Publications arising from IFIP events vary. The papers presented at the IFIP World 
Computer Congress and at open conferences are published as conference proceedings, while 
the results of the working conferences are often published as collections of selected and 
edited papers. 

Any national society whose primary activity is in information may apply to become a full 
member of IFIP, although full membership is restricted to one society per country. Full 
members are entitled to vote at the annual General Assembly, National societies preferring 
a less committed involvement may apply for associate or corresponding membership. 
Associate members enjoy the same benefits as full members, but without voting rights. 
Corresponding members are not represented in IFIP bodies. Affiliated membership is open 
to non-national societies, and individual and honorary membership schemes are also offered. 




Control of Distributed 
Parameter and 
Stochastic Systems 



Proceedings of the /F/P WG 7.2 
International Conference , June 19-22, 1998 
Hangzhou, China 



Edited by 

Shuping Chen 

Zhejiang University 

Xunjing Li 

Fudan University 

Jiongmin Yong 

Fudan University 

Xun Yu Zhou 

The Chinese University of Hong Kong 




SPRINGER SCIENCE+BUSINESS MEDIA, LLC 




Library of Congress Cataloging-in-Publication Data 

A C.I.P. Catalogue record for this book is available 
from the Library of Congress. 



ISBN 978-1-4757-4868-0 ISBN 978-0-387-35359-3 (eBook) 

DOI 10.1007/978-0-387-35359-3 

Copyright ® 1999 by Springer Science+Business Media New York 
Originally published by Kluwer Academic Publishers in 1999 



All rights reserved. No part of this publication may be reproduced, stored in a 
retrieval system or transmitted in any form or by any means, mechanical, photo- 
copying, recording, or otherwise, without the prior written permission of the 
publisher, Springer Science+Business Media, LLC 



Printed on acid-free paper . 




Table of Contents 



Preface vii 

Participants of the Conference ix 

Part I. Distributed Parameter Systems 

G. Avalos and I. Lasiecka: Exact- Approximate Boundary Controllability 

of Thermoelastic Systems under Free Boundary Conditions 3 

M. Bergounioux and F. Troltzsch: A Linear Parabolic Boundary Control 

Problem with Mixed Control-State Constraint 13 

M. Delfour: Membrane Shell Equation: Characterization of the Space of 
Solutions 21 

R. H. Fabiano: Renorming for Elastic Systems with Structural Damping . . 31 

F. Fahroo and C. Wang: Stability and Approximation of an Acoustic- 

Structure Model 39 

S. Hansen and Z. Liu: Analyticity of Semigroup Associated with a 

Laminated Composite Beam 47 

S. Kang, T. B. Stauffer, and K. Hatfield: A Practical Estimation Technique 
for Spatial Distribution of Groundwater Contaminant 55 

J. Lagnese: Domain Decomposition in Optimal Control of Elliptic Systems 

on 2-d Networks 63 

I. Lasiecka, R. Triggiani, and P. Yao: An Observability Estimate in 

1/2 (Q) x i7 _1 (Q) for Second-Order Hyperbolic Equations with Variable 
Coefficients 71 

S. M. Lenhart, M. Liang, and V. Protopopescu: Identification Problem for 

a Wave Equation via Optimal Control 79 

X. Li: Optimal Control Theory: from Finite Dimensions to Infinite 

Dimensions 85 

K. Liu and Z. Liu: Boundary Stabilization of a Hybrid System 95 

K. Liu and D. L. Russell: New Meaning of Exact Controllability of Linear 

Systems in Hilbert Spaces 103 

B. Mordukhovich: Minimax Design of Constrained Parabolic Systems .... Ill 

T. Nambu: Stabilization of Linear Boundary Control Systems of Parabolic 

Type: An Algebraic Approach 119 

T. I. Seidman: A Distributed Bioremediation Problem with Modal 

Switching 127 

G. Wang: Optimal Control Problems Governed by an Elliptic Differential 

Equation with Critical Exponent 133 

M. Yamamoto: Reconstruction of Source Terms in Evolution Equations 

by Exact Controllability 143 

J. Ye: Necessary Optimality Conditions for Control of Strongly Monotone 

Variational Inequalities 153 

Z. H. Zhang and X. Xiang: Optimal Controls of a Class of Strongly 

Nonlinear Evolution Systems 161 




VI 



Part II. Stochastic Systems 

E. K. Boukas and H. Yang: Robust Stabilization of Nonlinear Systems 

with Markovian Jumping Parameters 173 

S. Chen and X. Yu: Linear Quadratic Optimal Control: from Deterministic 

to Stochastic Cases 181 

P. Codings and U. G. Haussmann: Optimal Portfolio Selection with 

Transaction Costs 189 

T. Duncan: Some Approaches to Ergodic and Adaptive Control of 

Stochastic Semilinear Systems 199 

Y. Fujita: A One-Dimensional Ratio Ergodic Control Problem 207 

M. James: Nonlinear i/oo Control: A Stochastic Perspective 215 

M. Kohlmann: Reflected Forward Backward Stochastic Differential 

Equations and Contingent Claims 223 

H. Kunita: Short Time Asymptotics of Random Heat Kernels 231 

J. Ma and T. Zajic: Rough Asymptotics of Forward-Backward Stochastic 

Differential Equations 239 

J. B. Moore, X. Y. Zhou, and A. E. B. Lim: On LQG Control of Linear 

Stochastic Systems with Control Dependent Noise 247 

H. Morimoto and Y. Fujita: Radial Symmetry of Classical Solutions for 

Bellman Equations in Ergodic Control 255 

S. Peng: Open Problems on Backward Stochastic Differential Equations . 265 

R. Situ: Comparison Theorem of Solutions to BSDE with Jumps, and 

Viscosity Solution to a Generalized HJB Equation 275 

S. Stojanovic: Multivariate Constrained Portfolio Rules: Derivation of 

Monge- Ampere Equations 283 

L. L. Xie and L. Guo: Limitations and Capabilities of Feedback for 

Controlling Uncertain Systems 291 

G. Yin and M. Kniazeva: Time-scale Separation and State Aggregation 

in Singularly Perturbed Switching Diffusions 299 

J. Yong: Stochastic Controls and FBSDEs 307 

Q. Zhang and G. Yin: Asymptotically Optimal Controls of Hybrid LQG 

Problems: Summary of Results 315 

X. Y. Zhou and D. Li: Explicit Efficient Frontier of a Continuous-Time 

Mean-Variance Portfolio Selection Problem 323 




Preface 



The Conference on Control of Distributed Parameter and Stochastic Systems 
was held in Hangzhou, China, June 19-22, 1998. It was devoted to the current 
issues on the controls of distributed parameter and stochastic systems. The 
objective of the conference was to provide a forum where leading researchers 
from around the world can converge to disseminate new ideas and discusses 
the latest as well as the future trends in the areas. The following subjects were 
covered: Adaptive Control, Controllability, Filtering, Identification, Manufac- 
turing Systems, Mathematical Finance/Insurance, Numerical Approximation, 
Optimal Control, Stabilization, and Stochastic Analysis. 

There were two parallel sessions, distributed parameter systems and stochas- 
tic systems, in the conference. Totally 52 invited speakers delivered their 40- 
minute lectures. Participants came from the following 9 countries: Australia, 
Canada, China, France, Germany, Japan, Korea, Spain, and USA. 

In the distributed parameter session, the following were the main topics: (1) 
optimal control for PDE systems (by J. Burns, E. Casas, H. Gao, J. Lagnese, 
S. M. Lenhart, X. Li, T. I. Seidman, F. Troltzsch, G. Wang, X. Xiang, and 
J. Ye). (2) controllability /stabiliz ability of PDE systems (by M. Delfour, S. 
Hansen, M. A. Horn, I. Lasiecka, K. Liu, Z. Liu, T. Nambu, R. Triggiani, and 
E. Zuazua). (3) numerical study/ approximation (by R. H. Fabiano, F. Fahroo, 
S. Kang, C. Wang, and J. Zhou). 

In the stochastic session, the following two topics were very strongly repre- 
sented: (1) backward stochastic differential equations (by M. Kohlmann, J. Ma, 
S. Peng, R. Situ, S. Tang, and J. Yong), and (2) mathematical finance/insurance 
(by T. Bielecki, N. El Karoui, U. Haussmann, S. Stojanovic, M. Taksar, and 
X. Y. Zhou). Other topics included: linear-quadratic controls (S. Chen, J. 
B. Moore, Q. Zhang, and X. Y. Zhou), risk-sensitive controls (M. James and 
H. Nagai), ergodic controls (Y. Fujita, H. Morimoto and H. Nagai), adaptive 
controls (L. Guo and T. Duncan), Zakai equations (H. Kunita), singular per- 
turbation (G. Yin), and robust stabilization (H. Yang). 

The conference was financially sponsored by the Education Ministry of China, 
the National Natural Science Foundation of China, Zhejiang University and 
Fudan University. We also acknowledge the International Federation for Infor- 
mation and Processing (IFIP). 

This volume contains 39 written version of the talk presented at the confer- 
ence. 
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EXACT-APPROXIMATE BOUNDARY 
CONTROLLABILITY OF 
THERMOELASTIC SYSTEMS UNDER 
FREE BOUNDARY CONDITIONS 
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Abstract: Controllability properties of a partial differential equation (PDE) 

model describing a thermoelastic plate are studied. The PDE is comprised of a 
Kirchoff plate equation coupled to a heat equation on a bounded domain, with 
the coupling taking place on the interior and boundary of the domain. The cou- 
pling in this PDE is parameterized by a > 0. Control is exerted through the 
(two) free boundary conditions of the plate equation, and through the Robin 
boundary condition of the temperature. These controls have the physical in- 
terpretation, respectively, of inserted forces and moments, and prescribed tem- 
perature, all of which act on the edges of the plate. The main result here is 
that under such boundary control, and with initial data in the basic space of 
wellposedness, one can simultaneously control the displacement of the plate ex- 
actly . , and the temperature approximately. Moreover, the thermal control may 
be taken to be arbitrarily smooth in time and space, and the thermal control 
region may be any nonempty subset of the boundary. This controllability holds 
for arbitrary values of the coupling parameter a. 
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1 INTRODUCTION 



Statement of the Problem 



Let Q be a bounded open subset of M 2 with its sufficiently smooth boundary 
denoted as T. The boundary will be decomposed as T = To U IT, with both 
To and T i being open and nonempty, and further satisfying To H IT = 0. In 
addition, let V 2 be any open and nonempty subset of V. With this geometry, 
we shall consider here the following thermoelastic system on finite time (0 , T): 



o„(0,T)xfi; 

f 39 t — rjAd + a6 — a Ac o t = 0 v ; 



r\ 

w= ^ = 0 on (0,T) x r 0 ; 



Aw + (1 — + ad = ui 

< 9 Aw „ s dB 2 u du tt dd on( 0 ,r)xr i; 



( 1 . 1 ) 



dB +\d-{ u 3 °n (0, T) x r 2 

ai/ + 1 o on(o,T)xr\r 2 - u ’ 



cj(t = 0) = w 0 , = 0) = cji, 0(* = 0) = 6 0 on Q. 



Here, a, /?, rj and <r are positive parameters. The positive constant 7 is propor- 
tional to the thickness of the plate and assumed to be small with 0 < 7 < M. 

d 2 (jj 9 d 2 uj 9 d 2 uj 

The boundary operators Bi are given by B = 2z^i ^2 -7^ — ^ ^ — ^9 

axay ay 2 



d 2 uj 
dy 2 

the familiar Poisson’s ratio, and 1/ = [^1,1^2] denotes the outward unit normal 
to the boundary. Here we shall also make the following geometric assumption 
on the (uncontrolled) portion of the boundary To: 




and B 2 u = (v 2 - ^2) + v\vi 



3 {zq, yo} 6 M 2 such that /i(:c, y) • 1/ < 0 on To, (1.2) 



where h(x t y) = [x - x Q} y - y 0 ]. 

The PDE model (1.1), with boundary functions u\ = u 2 = 0, and U3 = 0, 
mathematically describes an uncontrolled Kirchoff plate subjected to a ther- 
mal damping, with the displacement of the plate represented by the function 
u(t 1 x ) y) ) and the temperature given by the function 9 (t } x,y) (see [8]). The 
given control variables ui(t,x) and u 2 (t,x) are defined on the portion of the 
boundary (0 ,T) x r 1; the control 1*3 (t,a:) is defined on (0 , T) x T 2 . 

With the denotation 



#P o (fi) = < tu E #*(0) : 



dvi 



— 0 for j — 0, ..., k — 1 
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we will throughout take the initial data #o] to be in #p 0 (ft) x #p o (ft) x 

L 2 (ft). In this paper, we will study controllability properties of solutions of 
(1.1) under the influence of boundary control functions in preassigned spaces. 
In particular we intend to address, on the finite time interval [0, T], the following 
problem of exact-approximate controllability with respect to the basic energy 
space JY 2 o (ft) x #p o (Q) x T 2 (ft) (this term being coined in [5]): For given data 
[u>o , , 0 Q ] (initial) and [loq , ujJ , ] (terminal) in H^ 0 (ft) x H^ 0 (ft) x L 2 (ft) and 

arbitrary e > 0, we seek a suitable control triple [ 14 , 1 * 2 , 1 * 3 ] £ L 2 (0,T] L 2 (Ti) x 
H-'iTx)) x C r (E 2| T) (where r > 0 and E i|T = (0,T) x r z , i = 0,1,2) such 
that the corresponding solution [uj^u} t) 9\ of (1.1) satisfies the steering property 
at terminal time T 

[w(T)MT)] = K,^f] i and ||0(T) - || ia(n) < (1-3) 

In regards to the literature on this particular problem, the most relevant 
work is that of J. Lagnese in [9]. Therein, Lagnese shows that if the coupling 
parameter a is small enough and the boundary T is “star-shaped”, then the 
boundary-controlled system (1.1) is (partially) exactly controllable with re- 
spect to the displacement u. Also in [14], a boundary-controlled thermoelastic 
wave equation is studied, with a coupling parameter a likewise present therein, 
and a result of partial exact controllability (again for the displacement) for this 
PDE is cited. This controllability result is quoted to be valid for all sizes of a; 
however in [15], the author has acknowledged a flaw in the controllability proof, 
the correction of which necessitates a smallness criterion on a, akin to the situ- 
ation in [9] . The chief contribution of the present paper is to remove restrictions 
on the size of the coupling parameter (see Theorem 3 below) , at the expense 
of adding the arbitrarily smooth boundary control us in the thermal compo- 
nent. For a 1-D version of (1.1), S. Hansen and B. Zhang in [6], via a moment 
problem approach, show the system’s exact null controllability with boundary 
control in either the plate or thermal component. Other controllability results 
for the thermoelastic system which do not assume any “smallness” condition on 
the coupling parameters deal with distributed /internal controls. Such include 
that in [5], in which interior control is placed in the Kirchoff plate component 
subject to clamped boundary conditions; with such control, one obtains exact 
controllability for the displacement u>, and approximate controllability for the 
temperature 6. Alternatively in [3], interior control is placed in the heat equa- 
tion of (1.1) so as to obtain exact controllability for both components u and 6. 
In addition, the work in [13] deals with obtaining a result of null controllability 
for both components of a coupled wave and heat equation, in the case that 
interior control is inserted in the wave component only. 

So again, the main contribution and novelty of this paper is that we consider 
boundary controls acting via the higher order free boundary conditions, and we 
do not assume any size restriction on the coupling parameter a. Moreover, 
we do not impose any geometric “star-shaped” conditions on the controlled 
portion of the geometry. 
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It should be noted that the particular type of boundary conditions imposed 
on the mechanical variables greatly affects the analysis of the problem, even 
in the case of internal control. Indeed, in the case of all boundary conditions, 
save for the free case , it is known that the thermoelastic plate semigroup can be 
decomposed into a damped Kirchoff plate semigroup and a compact perturba- 
tion (see [5] and [12]). Since controllability estimates are invariant with respect 
to compact perturbations (at least in the case of approximately controllable 
systems, which we are dealing with here), the aforesaid decomposition, valid 
for the case of lower order boundary conditions, reduces the problem of exact 
controllability for the mechanical variable to that of uncoupled Kirchoff plates. 
Thus, the case of lower order boundary conditions allows a reduction of the 
coupled problem into one which has been much studied in the past. This strat- 
egy, while successfully employed in the case of clamped or hinged boundary 
conditions (see [5]), is not applicable here. Indeed, in our present case of free 
boundary conditions, there is no decomposition with a compact part, as in the 
lower order case (see [ 12 ]); moreover, the controllability operator correspond- 
ing to the given boundary controls is not bounded on the natural energy space. 
This latter complication is due to the fact that the Lopatinski conditions are 
not satisfied for the Kirchoff model under free boundary conditions. 

The strategy adopted in this paper consists of the following steps: Initially, a 
suitable transformation of variables is made and applied to the equation ( 1 . 1 ); 
subsequently, a multiplier method is invoked with respect to the transformed 
equation. The mulitiplers employed here are the differential multipliers used in 
the study of exact controllability for the Kirchoff plate model, together with the 
nonlocal (^DO) multipliers used in the study of thermoelastic plates in [ 1 ] and 
[2]. This multiplier method allows the attainment of preliminary estimates for 
the energy of the system. However, these estimates are “polluted” by certain 
boundary terms which are not majorized by the energy. To cope with these, we 
use the sharp trace estimates established in [11] for Kirchoff plates. The use of 
this PDE result introduces lower order terms into the energy estimate, which 
are eventually eliminated with the help of a new unique continuation result in 
[7] . It is only at the level of invoking this uniqueness result that the thermal 
control U 3 on T 2 must be introduced. 

We post our main result here on controllability. 

Theorem 1 Let the assumption (1.2) stand. There is then a T* > 0 so that 
for T > T* the following controllability property holds true: For given initial 
data [w o ,wi, 0 o] an d terminal data [(Jq , uif , 9q] in the space H^ q (Q) xi7p o (Q) x 
L 2 (Q), and arbitrary e > 0 , one can find control functions [u*, 11 $, U 3 ] E 
L 2 (0, T; T 2 (ri) x H~ 1 ( Ti)) x C r (£ 2 ,T) (where given r > o) such that the cor- 
responding solution [cj*,cj*,0*] to (1.1) satisfies (1.3) at terminal time T. 

Remark 2 The presence of the control in Theorem 1 is owing solely to 
the need to invoke the aforementioned uniqueness result of Isakov in the proof 
below; it plays no part whatsoever in obtaining the preliminary (lower order 
term-tainted) estimate on the energy. Consequently, we have the freedom to 
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prescribe the thermal control region to be as small as we wish , and the control 
113 to be as smooth in time and space as desired . 

2 PROOF OF THEOREM 1 

A preponderant portion of the proof of Theorem 1 is wrapped up in showing 
the following result of exact controllability for the displacement only: 



Theorem 3 With the coupling parameter a in (1.1) being arbitrary and the as- 
sumption (1.2) in place , there is then aT* > 0 so that for T > T* , the following 
property holds true: For all initial data [^ 0,^1 > 0o] € H^ q (Q) x (0) x L 2 (Q) 
and terminal data E H p (Q) x H p ( 0 ), there exists [111,1x2)^3] E 

L 2 ( 0 , T; L 2 (Ti) x H~ 1 ( Ti)) x H s (Yi2 ) t)> where arbitrary s > 0 , such that the 
corresponding solution [w,tj t ,0] to (1.1) satisfies [w(T), u t (T)] = . 



Indeed, if Theorem 3 is shown to be true, then using the miminal norm 
steering control (see Appendix B of [10]), one can, in a straightforward fash- 
ion, construct a control [u*, txj, ixj] such that the corresponding trajectory 
[u*,l <;*,$*] has the desired reachability property (1.3). (See [4] for the pre- 
cise details). Accordingly, the sequel is devoted to showing the validity of 
Theorem 3. 

The theme of the proof of Theorem 3 is a classical one. Denoting the 
control space U s = T 2 (Ti) x x H S (T 2 ), one defines the operator Ct : 

D(Ct) C U — > H p (fi) x f7p o (Q) to be that which takes the terminal control 
to the terminal state; i.e., Ct [ 1 x 1 , 1 x 2 , 1 x 3 ] = M T),u t (T)], where [u(T),u t (T)] 
is the plate component of the solution to (1.1) at time t — T. As it is defined, 
Ct is a closed, unbounded operator, with its domain being densely defined. 
By a principle of functional analysis then (see e.g., [16]), to prove Theorem 
3, which is essentially a statement of the surjectivity of it is enough to 
establish the PDE inequality 




2 

L 2 (r 1 ) 



^+IMI [ h s { z 2 , t )]' - Ct ||[0O)^i]|| 



2 



( 2 . 1 ) 



where [0, <j> t) ip\ is the solution to the following backwards system, corresponding 
to terminal data [<j> 0 , <j) 1 ] E D(Ct ) • 
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f <f) tt - 7 A <p t t + A 2 <p + a Aip = 0 
\ /3ip t + rjAip - <np - a A <p t = 0 



on (0, T) x 



dp 

p = — = 0 on £ 0 ,t; 



A</> + (1 - /i)5i0 + aip — 0 



9A0 dB 2 <p 

-sr + (1 -' ,) 



d<j>tt dip _ on Si,t; 



_ 7 _l + a ^ = 0 

dr 7 fli/ + <9z/ 



dip 

du 



+ = 0 on (0, T) x T, A > 0; 



fe(T),^T),tf(T)] = [^,^,0] 



( 2 . 2 ) 



The proof of Theorem 3 is then based upon showing the inequality (2.1), 

dpt I 



at least for T > 0 large enough, where 



dv 



Pt\r x 5 ^1 



r 2 



are traces of the 



solution [<p ) <p t , ip] to the backwards system (2.2). Because of space constraints, 
we give here only a broad sketch of the proof of Theorem 3; the full particulars 
are provided in [4]. 

Step 1. We start by making the substitution 



<p(t) = e ^(p{t); and ip(t) = e ^ ip(t ), (2.3) 



2 

where parameter £ = 5—. This particular choice of parameter allows the PDE 
(2.2) to be transformed into the following plate equation whose forcing function 
is comprised in part of the high order term ip t : 



(ptt ~ lA(p tt -f A 2 (p = c Q ip + ciipt + c 2 p + c 3 <p t + c 4 A<p on (0, T) x Q; 



- dp 

<P = = 0 on E 0)T ; 






( A<p + (1 — ii)B\<p — — aip 
{ dAj , ,, dB 2 p d ( 

\ 0 , +(1 "> Sr =1 eA 


£ 2 P + 2£(pt -f Ptt'j 


di> 

~ a !& 



= [e-^ o ,-£e-^ o+ e-^ 1)0 ] , 



on Ei^; 



(2.4) 



where the c; are constants depending on the physical parameters. 

Step 2. At this point we invoke a multiplier method with respect to the 
uncoupled equation in (2.4), using two multipliers which are, respectively, stan- 
dard (as in [9]) and nonstandard (as in [1] and [2]). To wit, we multiply this 
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equation by <j > t , and subsequently integrate in time and space so as to arrive at 






rT L, 



Iff* (n) 



= ~ a l 



r / d* t 

^>- 5 — 

ov 



- I 7^ + + I dt 

Js \ " /L=(r.) 

+ / (coV’ + c 2 ^+ c^ t + c 4 A^+ ciV’t, ) dt. (2.5) 

J s v y i 2 (ri) 

Moreover, letting Ad denote the Laplacian operator with Dirichlet boundary 
conditions, and Af, 1 £ £(L 2 (ft)) its corresponding (smoothing) inverse, we then 
multiply the PDE in (2.4) by -^-Af^xj), and thereafter integrate in time and 
space. Adding the resultant expression to (2.5), and subsequently majorizing 
the sum, we obtain the following: 



Lemma 4 The solution to (2-4) satisfies the following relation for 

all s and r £ [0, T]: 






+ N*)L 






nT 2 

c / V^t dt+ l.o.t. , 

Jo ^(ro ' ' 



where l.o.t yip } <t>, <l>tj denotes, as usual, “ lower order terms ” (below the en- 
ergy level) ofip, <f), and < j> t . 

Step 3. Taking a radial vector field h G M 2 which meets the requirement 
in (1.2), one can derive the following inequality which is an analogue to that 
demonstrated in [9]. In deriving this estimate the trace result in [11] is critically 
invoked. 

Lemma 5 For all €q E (0,T), the solution to (2.4) satisfies 



rT-e 0 ^ 2 xv 2 

/ 0 + <t*t 

Je o [ *? 0 (n) *J 0 («) 



< C T / V& 



fa { s i , 






where s \ — T — c 0 , and s 2 = eo- 
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Step 4. A standard energy argument shows that the heat component ip 
satisfies the estimate 




This inequality, in combination with Lemmas 4 and 5, eventually give (again 
the full details are in [4]), 



Lemma 6 ForT > 0 large enough, the solution of (2.4) satisfies the 

following estimate: 



< 






| 2 




2 




2 




T 


1 




Pp 






1 *?„(«) 




if 1 (ft) 



+ 



kT),MT) 



2 





i2(ri) *+ 1-o.t. (£&,£)• 



( 2 . 6 ) 



Step 5. Note that in the inequality (2.6), there is no boundary trace term 
ip Ip 2 , reflecting the contribution of the control u$\ the observability estimate 
(2.6) is independent of thermal control. However, to remove the corrupting 
lower order terms in this estimate so as to have the desired inequality (2.1), 
the thermal control now comes directly into play. Indeed, a compactness- 
uniqueness argument is now to be employed, with this argument making critical 
use of the Holmgren’s-type uniqueness result derived in [7] for overdetermined 
(in both the mechanical and thermal variables) thermoelastic systems. The cor- 
rect use of this uniqueness theorem necessitates the appearance of the thermal 
control U3 . With such control in place, we then have 



Lemma 7 ForT > 0 large enough, the existence of the inequality (2.6) implies 
that there exists a Ct such that the following estimate holds true : 



I, 



/-V — - 


~\ ( r T 




2 




\4>AtA 


’) s Ct u 




dt T 

L2 ( r 0 





[tf s (S 2)T )]' 



(2.7) 



The inequalities (2.6) and (2.7), and the transformation <p{t) — e ^(p(t) 
and ip(t) = e~^ip(t) give the desired inequality (2.1), thereby completing the 
proof of Theorem 3. By the remarks made at the beginning of this section, 
with this partial exact controllability result in hand, the exact- approximate 
controllability statement Theorem 1 follows. 
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Abstract: 

A simple class of linear optimal control problems for parabolic equations 
with mixed control-state inequality constraints is investigated. The constraints 
are formulated pointwise in L°°. It is shown how to obtain associated Lagrange 
multipliers in L p - spaces. 

1 INTRODUCTION 

We discuss the following linear optimal boundary control problem for the heat 
equation: 

max / an y(T) dx + olq y dxdt + / a^y dadt- h / a u u dadt 

Jn Jq Jt. */£ 

subject to the state equation 

y t - Ay +dy = 0 inQ, 

d„y + by — u on E , (1.1) 

y(0) = 0 in £1 , 

and to the mixed control-state constraints 

u(x,t) < c(x,t) + y(x,t), 
u(x,t) > 0, 

which are required a.e. on E. The heat equation is defined in a bounded 

domain C IR^, N > 2, with sufficiently smooth boundary T. In the fixed 
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time interval (0 , T) we put Q := Q x (0 , T), and E := T x (0,T). Moreover, 
functions G L°°(0), d, olq G L°°(Q), and c, 6, G L°°(E) are given. 

The control function u is assumed to be bounded and measurable. In this way, 
the feasible set of the control problem belongs to L°° (E), and the mixed control- 
state constraint u < c + y must be regarded in the same space. Therefore, one 
might expect that an associated Lagrange multiplier has to be found in L°° (E)* . 
In contrast to this, we shall verify the existence of at least one multiplier in 
L°°(E). For linear programming problems in L p -spaces with constraints of 
bottleneck type this surprising fact is known since long time, see for instance 
[1], [2], and the references cited therein. In this short note, we extend these early 
ideas to the parabolic boundary control problem defined above. Asscociated 
distributed control problems have been discussed extensively in our recent paper 
[3] . In the forthcoming paper [4] , a class of nonlinear parabolic control problems 
with pointwise mixed control-state inequality constraints will be discussed on 
using these results. 

We assume that T is so smooth that a Green’s function G — G(x ) ^ ) t) ) 
G : Q 2 x IR.f IR exists, which satisfies the known estimate 

|G(g,£,<)| < kit~* exp(— k 2 (1-2) 

with positive real numbers k\ ) Then the (weak) solution of (1.1) is given 

by 

t 

y{x,t) = J J G(x,£,t- s)u{Z,s)d<r(Z)ds, (1.3) 

0 s 

where da denotes the surface measure on V. 

2 COMPARISON PRINCIPLES FOR AN INTEGRAL EQUATION 

First we discuss the integral equation 

u(x,t) - f{x,t) + j J G(x,Z,t-s) u(£,s)da(£)ds (2.1) 

for given / G L°°(E). Introduce the integral operator K: 

( Ku)(x,t ) — a G(x,£,t - s) u(£,s) da(£)ds. (2.2) 

I\ is continuous in L°°(E), also from L p ( E) to L°°(E) for p > N + 1, see [5], 
p. 138, Lemma 5.6.6. K is also continuous from L P (E) to C(Q). In this way, 
2 /(-,T) is well defined in C(Cl). Endow the space L°°(E) with the equivalent 
norm ||u||x, 

||u||a = vraimax|e _At u(x,^)| , (2.3) 

(a7,t)(E£ 

where A > 0. Then it is an easy and standard exercise to show the following 
result: 
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Lemma 2.1 K is a contraction in T°°(E), provided that A > A 0 is sufficiently 
large. 

Corollary 2.1 For all f E L°°( E), the equation (2.1) has a unique solution 
u E L°°(E). The mapping f u is continuous in L°° (1 1). IfG(x,£,t) > 0 and 
f(x,t ) > 0 a.e. in E, then u(x,t) > 0 a.e. in E. 

Proof: We have 

oo 

« = (/- K)~ 1 f = (J2 K n )f (2.4) 

n=0 

by well known results on Neumann series. Moreover, 1 1 (/ — isT) ~ 1 1 | a < 1/(1 — 
||/i ||a), where we have used for convenience the symbol || • \\\ to denote also the 
norm of the operator K induced by || • 1 1 x • The first part of the lemma follows 
from the equivalence of the norms || • ||a and || • ||l°°(e)- If G is a nonnegative 
function, then K is a nonnegative operator, that is / > 0 => K f > 0. The 
second result follows immediately from this and (2.4). ■ 

If /i , /2 are two bounded and measurable right hand sides for (2.1) and G 
is nonnegative, then u\ > U 2 holds for the associated solutions. Now we verify 
these facts for the space L p (Yf). To see this, regard K as an operator from 
L p ( E) to L°°(E). This smoothing property of K is basic for the next result. 

Corollary 2.2 (Comparison principle) Suppose p> N + 1. Then the inte- 
gral equation (2.1) has for each f E L p { E) a unique solution u E T P (E), and the 
mapping f i-» u is continuous in L P (E). IfG(x,£,t) > 0, Ui E L p ( E), i = 1,2, 
are solutions of (2.1) associated to fi E L p ( E), i — 1,2, and fi(x } t) > f 2 (x,t) 
holds a.e. in E, then ui{x ) t) > U 2 (x,t) a.e. m E. 

Proof: Put v := u — /, then (2.1) reads v — Kf + Kv. By the smoothing 
property of K, the right hand side I< f of this transformed equation is bounded 
and measurable. The last corollary implies that this equation admits a unique 
solution v E L°°(E) depending continuously on K f and hence on /. Clearly, 
u = / -f v is a solution of (2.1) in L p ( E). The uniqueness of u in L p ( E) is 
and easy consequence, since the difference of two solutions solves the equation 
with right hand side zero, which belongs to L°°(E). The comparison part of 
Corollary 2.2 follows from the arguments after Corollary 2.1. ■ 

In a dual problem, which is defined later, the (formal) adjoint integral oper- 
ator A' t , 

(K T n){x,t) = ( jG(Z,x,s-t)n(£,s)dcr(Z)ds (2.5) 

appears. K J has the same properties as K. In particular, for A > \ 0 it is a 
contraction in the norm || • 1 1 a • 

Lemma 2.2 For every function a E T°°(E), the equation 

p{x ) t) — max{ 0, a(x,t) -f a G(£, x,s-t) s) da{£)ds} (2.6) 




16 



has exactly one solution p E Z/°°(E). 

The proof is an application of the contraction mapping principle. It can be 
applied, since the operator (II z)(x,t) = max{ 0, z(x, t) } is Lipschitz continuous 
with Lipschitz constant one in L°°(E), and K J is a contraction. 



3 PRIMAL AND DUAL PROBLEM 

By inserting the integral representation (1.3) of the state y in the objective 
functional of the control problem, we get after changing the order of integration 

/ any(T)dx-\- aQydxdt-\- y dadt-\- a u u dadt — (—a(x ) t))u(x ) t)dadt 

Jfl J e J s is J 



with a certain function a E L°°(E) (a concrete expression for a is given in 
section 5). In this way, the optimal control problem becomes a continuous 
linear programming problem with constraints of bottleneck type. This is our 
Primal Problem 



max 



L a 



dadt 



CP) { 



u(x,t) < [c 4* Ku](x,t) a.e. in E , 
u(x,t) > 0 a.e. in E . 



Note that (V) is defined in the space L°°(E). From now on, we assume that 
G(x, £, t — s) > 0 and p > N -f 1. We define u c by the equation 



u c — c + Ku c 



(3-1) 



According to Corollary 2.1, u c is bounded and measurable. 



Theorem 3.1 (V) has a solution u if and only if u c > 0. 



Proof: If u c > 0, then the feasible set for (V) is non-empty. All feasible elements 
u E L p ( E) satisfy u < u c by the comparison principle of Corollary 2.2, and we 
have u c E L°°(E). Moreover, u > 0 follows from the constraints. This implies 
the L°°-boundedness of u by ||i/c||l oo (E)- So the feasible set is bounded, closed, 
and convex in all (reflexive) L p -spaces for 1 -f N < p < oo. The existence of 
u E L p ( E) is an immediate conclusion. Of course u belongs to L°°(E). On the 
other hand, if u c is not greater or equal than 0, then the feasible set is empty 
by the comparison principle. ■ 

To introduce the dual problem to (V), we extend the feasible set of (V) 
from L°°(E) to L p ( E) (p > N + 1). Any feasible solution u of (V) satisfies 
0 < u < u c , where u c is defined in (3.1). Therefore, all feasible solutions of 
L P (E) belong automatically to L°°( E), and the feasible set is not influenced by 
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this extension to L P (E). From standard techniques to establish dual problems 
(see our discussion in section 6), in L 9 (E) we get the following Dual Problem : 

min J c p dadt 

p{x,t) > [a + K r p\(x,t) a.e. in E , 

p(x, £) > 0 a.e. in E , 

where p E L 9 (E), — I — = 1, and K J is the adjoint integral operator defined 
P Q 

by (2.5). The kernel of K J satisfies the estimate (1.2), hence K r : L p ( E) — > 
L°°(E) holds for p > N + 1, too. However, this is not true from L g (E) to 
L°°(E). On the other hand, K r represents the adjoint operator of K : L P (E) — ► 
L°°( E) C L p (Yj). Therefore, K J is continuous in L q (E) ~ L P (E)*, too. 

In view of (1.2) and Lemma 2.1, K T is a contraction in L°°(E), hence the 
equation 

n = p + K T II (3.2) 

has a unique solution [i E L°°(E) for each /3 E L°°(E). Moreover, we have 
uniqueness for (3.2) in L q ( E). This can be shown by duality, since by Corollary 
(2.2) the range of v — Kv is L p ( E). The next theorem contains the basic idea 
of this paper. 

Theorem 3.2 If c > 0, then (V) has at least one bounded and measurable 
optimal solution p. 

Proof: We know by Lemma 2.2 that equation (2.6) 

p(x,t) = max{ 0, a(x,t) + (A" T u)(£,t)} 

has in L°°(E) exactly one solution fi. Let p E L q ( E) be any other feasible 
element for ( V ), which is different from p. Then p > a + K T p a.e. in E and, 
of course, p > 0. Next, we construct a sequence p\ > p 2 > • • • as follows : 
pi = p E L q (Yf) and p 2 = max{ 0, a + K J pi } . A simple discussion yields 
P2 < Pi a.e. in E. Then by positivity, K J p 2 < K T pi, and we get 

p 2 > a + AT T /ii > a -|- K T p 2 a.e. on E . 

So /^2 is feasible and /i 2 < Pi on E. Repeating this process, one constructs 
a non-increasing feasible sequence {p n } which has to be pointwise convergent 
towards some p > 0, that is 

lim p n (x,t) — p(x,t) a.e. on E . 

n-*+oo 

An application of the Lebesgue dominated convergence theorem yields that p n 
tends to p in L 9 (E). Passing to the limit in p n — max{ 0, a + K J p n -i } gives 
in L g (E) 

p — max{ 0, a -|- K T p } . 





18 



Finally, by simple arguments, we are able to conclude jl G L°°(£). By unique- 
ness in L°°(£) we have jl = jl] moreover p > jl. Therefore p > jl holds for all 
feasible solutions and, since c > 0, 



c(x,t) p(x,t) dadt > / c(x,t) p(x,t) dadt. 



4 THE DUALITY RELATION 

The discussion of the duality between (V) and (V) is not yet complete. We have 
only shown that under certain assumptions the dual problem admits a solution 
jl. To make sure that jl is a Lagrange multiplier associated to a solution of 
(‘ V ), we need additionally the strong duality relation , that is the equality of 
primal and dual optimal value. To this aim, we briefly sketch some main ideas 
of duality for linear programs. 

Let X = L p (£) with its natural partial ordering >, and define A := I — K 
having the adjoint operator A! . By (• , •) we denote the pairing between X and 
its dual space X' — L q ( £). Then the primal problem is 





max 


(a,x) 


(V) 


A x 


< c 




X 


> 0 . 



On using the Lagrange function C{x,p) := ( a,x ) + (p,c — Ax), the primal 
problem can be written in sup-inf form. Reversing the order of supremum and 
infimum we arrive at the dual problem 

min (//,c) 

(V) A! p >a 

/ i >0. 

Let x be optimal for (V). Then jl is an associated Lagrange multiplier if and 
only if the pair (x, jl) is a saddle point of C , i.e. C(x , jl) < C(x , jl) < C(x , p) for 
all x > 0, p > 0. A necessary and sufficient condition for [x,jl) to be a saddle 
point is that x solves (' V ), jl solves (‘ D ), and the strong duality relation v = v f 
holds true. Let us define by P(c) and D(a) the feasible sets of (V) and (V), 
respectively, i.e. 

P(c) = {x G X\x > 0, A x < c}, D{a) = {p E X'\ p > 0 , A' p> a). 

It is easy to verify that the weak duality relation 




holds always true. To show the strong duality relation, we need the convex 
cone 

K(A) := {(a, d) G IR x X\3x > 0 : Ax < d, (a, x) > a}. 
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Theorem 4.1 If K (A) is closed and (V) admits at least one solution , then the 
strong duality relation holds, that is 

ma x(a,x)= inf (u,c). n o\ 

xeP( c y ' t*eD{\y 1 (4-2) 



This is a standard result of linear programming theory. We should mention 
that (4.2) implies the complementary slackness conditions 

{A x - c, p) = {A' p- a,x) — 0. (4.3) 

It turns out that the assumptions of the theorem are fulfilled for our problem: 

Lemma 4.1 I< (A) is closed for ( V ). 

We know from the proof of the Theorem 3.1 that the norm of u E P(c) is 
bounded by the norm of u c , which is bounded by Corollary 2.2. Therefore, this 
Lemma is easy to prove. 

5 NECESSARY OPTIMALITY CONDITIONS 

Finally, we prove that the solution p of (V) is a Lagrange multiplier for the 
state-constraint of the parabolic control problem. By definition, p is an associ- 
ated Lagrange multiplier , if u , together with p and an adjoint state p , satisfies 
the first order necessary optimality conditions : They consist of the adjoint 
equation 



-pt-Ap+dp = 


ag 


in Q 




d„p +bp = 


as 


on E , 


(5.1) 


p(T) = 


&n 


in Q , 





the variational inequality 

J {ot u 4- P + P){u - u) dcrdt > 0 Vu > 0, (5.2) 

E 

the complementary slackness condition 

[u — c — y) p — 0 a.e. in E, (5.3) 

and the nonnegativity condition p(x,t) > 0, which must be satisfied a.e. on E. 

To verify these conditions, we introduce some auxiliary functions: Put an = 
—an, as = —as, ag = — ag, a u — — a u . Moreover, we define and <p by 

-^t ~ + = ag —<p t — A(p + d<p — 0 

d^ + b 1 # — as d y <p + b<p — p 

*{T) = an <p{T) = 0 . 




20 



The function ip plays the role of the dual state. It holds K r p = <p, hence 
the constraints of (V) admit the form p > a + <p. Moreover, (4.3) yields the 
complementary slackness condition for the dual problem, 

(— a — ip -f p) u — 0 a.e. on E. (5.4) 

was defined to satisfy a = ^A-a u = ^f — a u . 

It is clear that we have to put p := — (\£ -{-</?). Then p solves (5.1). Moreover, 
the complementary slackness condition (5.3) is satisfied by (4.3). Finally, the 
variational inequality follows from 

f(a u + P + - u) = f(-a u + p + p)(u - u) = f (^ - a + P + p){u - u) 

E s E 

= f (-a - <p + p){u - u) = / (-a — p -f p)u - f (-a - <p + Jl)u > 0. 

E E £ 

In the last estimate, the relations — a — <p + p > 0 and (5.4) were used. The 
optimality conditions are verified. Altogether, we have found our final result: 

Theorem 5.1 The dual problem (D) has at least one bounded and measurable 
solution p. Let u be optimal for the linear parabolic boundary control prob- 
lem. Then p is a Lagrange multiplier associated with the mixed control-state 
constraint u < c + y. 
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Abstract: The existence of solution to membrane shell equation is studied 

in a bounded open connected domain w (Lipschitzian when w has a boundary 
7) in a C 1 ' 1 midsurface for homogeneous Neumann boundary conditions or 
homogeneous Dirichlet boundary conditions on a part 70 of 7. It is proved that 
its tangential part is solution of the reduced membrane shell equation in H 1 (w) N 
(resp. Hy 0 (u>) N ) unique up to an element of a finite dimensional subspace, while 
its normal component belongs to a weighed L 2 (u>) space by the pointwise norm 
of the second fundamental form. It is also shown that the reduced equation is 
equivalent to the equation for the projection onto the linear subspace of vector 
functions whose linear change of metric tensor is orthogonal to the second 
fundamental form of the midsurface. 

1 INTRODUCTION 

In recent papers ([8, 9, 11, 12]) it was established that the polynomial P( 2, 1) 
model is both pertinent and basic in the theory of thin shells. It was shown in 
[8] that its solution converges to the solution of a coupled system of variational 
equations. For the plate and the bending dominated shell it yields (as the 
thickness 2 h goes to zero) the membrane shell equation and the asymptotic 
bending equation. 

The first variational equation of the asymptotic coupled system coincides 
with the variational equation characterizing the asymptotic P(0, 1) model. It 
was shown in [8] that this equation decomposes into two equations: a first 
equation containing the Love-Kirchhoff group of terms and a second equation 
which coincides with the classical membrane shell equation. The detailed cor- 
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respondence with the covariant form of the classical membrane shell equation 
is given in [10]. The decomposition is achieved by variable elimination which 
results in the introduction of an effective compliance C e p associated with the 
initial three dimensional compliance C. 

In this paper the membrane shell equation is studied in a bounded open 
connected domain lj (Lipschitzian when lj has a boundary 7) in a C 1 ’ 1 sub- 
manifold of codimension one for homogeneous Neumann boundary conditions or 
homogeneous Dirichlet boundary conditions on a part 70 of 7 when 70 has non- 
zero Hjsr- 2 Hausdorff measure. This paper is a companion paper to [8] where 
the spaces of solution E° and E p corresponding to the respective asymptotic 
P( 0, 1) model and the membrane shell equation are defined as completions of 
appropriate quotient spaces. It gives a complete characterization of the space 
E p without extra condition on the second fundamental form. Such a charac- 
terization is currently available for the plate and uniform strong elliptic shells 
in [13, 5, 6, 7]. It also shows that we can always associate with the vector 
functions of the space E p a class of tangential components which turns out to 
be solutions of the reduced membrane shell equation. This reduced equation 
is also connected with a projection onto a linear subspace of elements of E p 
whose linear change of metric tensor is orthogonal to the second fundamental 
form. Another consequence of the characterization of E p is the fact that in the 
asymptotic convergence of the solution of the P(2, 1) model we now know that 
the tangential component of the displacement of the midsurface strongly con- 
verges in H 1 ^) 1 * and the normal component in a weighed L 2 {lj) space by the 
pointwise norm of the second fundamental form. The characterization given in 
this paper and the one of E 01 given in [9] for the P( 2, 1) model sharpen the 
abstract results of [8] . 

For N = 3, this extends to arbitrary D 2 b the available existence of solutions 
obtained by [1, 3] for g° = 0, homogeneous Dirichlet boundary conditions on 
the whole boundary, the special constitutive law C~ l e — 2/i£ + Atr el and 
the uniform ellipticity of the 2-dimensional C 2 -midsurface cj. However in the 
case of uniform elliptic shells uniqueness does not so far follows directly in an 
obvious way. The first existence and uniqueness result seems to be due to [13] 
under relatively strong conditions. For a domain lj with a C 3 boundary 7 
in an analytic midsurface, the existence and uniqueness of solutions (£p, t)°) in 
Hq(u }) 3 x L 2 (lj) was established by [6, 7]. The conditions were relaxed by [1,3]: 
the midsurface is of class C 2 and the boundary 7 is Lipschitz for the existence 
(midsurface C 5 and the boundary 7 of class C 4 for existence and uniqueness) . 

Notation and Background Material. The inner product in R^ and the 
double inner product in ^(R^ ; R^) (space of N x N matrices or tensors) are 
denoted as 

xy = Eili x i Vi , A- B = Eili £f=i Aij Bij . 

* M denotes the transpose of of an arbitrary k x m matrix M. 

In this paper the submanifold T of codimension one is specified as the bound- 
ary of a subset fi of R^ . It is assumed that lj is a bounded open subset of T 
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and that T is of class C 1,1 in a neighbourhood of u. This is equivalent to say 
that the algebraic distance function 6 of Q, is C 1,1 in that neighborhood. Its 
gradient V6 coincides with the unit exterior normal n to T on uj and its Hessian 
matrix D 2 b to the second fundamental form. Finally P d = I — n*n will denote 
the orthogonal projection onto the tangent plane to to ([ n*n]ij = ninj). For 
a detailed account of the intrinsic differential calculus on a C 1, ^submanifold, 
the reader is referred to the now available lecture notes [10, 8]. Finally it will 
be convenient to introduce the following notation for the decompositions of an 
N x N matrix r into its tangential and normal parts along u 

r p d Ji f p r p ; T nn rn • n, [t]r = r p + (Pm) *n -f n*(Prn) + r nn n *n 
and the spaces of symmetrical matrices 

Sym N d = {r G £( R N - R N ):* r = r} 

Sym^ d = {r £ Sym N :rn = 0} Vr e Syrrijv, r p £ Sym^. 



2 MEMBRANE SHELL EQUATION 

It was shown in [8] that the membrane shell equation can be obtained by de- 
composition of the variational equation of the asymptotic P(l, 0) model which 
also yields the typical group of terms occurring in the Love-Kirchhoff condi- 
tion. It involves an effective compliance C e p which retains the properties of 
the three-dimensional compliance C. So for the purposes of this paper it is 
convenient to start with the following assumption on the effective compliance. 

Assumption 2.1 The effective compliance is a linear bijective and symmetri- 
cal transformation of Sym^ such that there exists a constant a > 0 for which 

Mr p £ Sym£, C~pT P --r p >a\\T P \\ 2 

The membrane shell variational equation is given by: for all v° £ i/ 1 (w) iv 

[ C^e£(v°)-e£(v°)dT = e p (v 0 ) (1.1) 

J u> 

where the right-hand side is specified by a linear functional t p . Associate with 
£p the space 

V d M {v£L 2 (u,) N :v r £H 1 (u,) N }cH d ^ L 2 (u,) n (1.2) 

and define E p as the completion of the quotient space V j ker e p with respect 
to the norm associated with the inner product 

[ e r( u )'‘ £ ri v ) dF. 



(1.3) 
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Similarly for homogeneous Dirichlet boundary conditions on a part 70 of 7, 
denote by E p q the completion of the quotient space 

V 70 /ker e p , V l0 = {# 6 L 2 (u) N :vr G H}, 0 (w)* } 

with respect to the norm generated by the scalar product (1.3). By Assump- 
tion 2.1 on C e p, the bilinear term in (1.1) is continuous and coercive in E p . 

Theorem 2.2 Let Assumption 2.1 on C e p be verified. 

(i) Given t p E {E p ) f , that is there exists c > 0 such that for all v° E H 1 (oj) n 

rK)|<c|| £ f(,°)|| L2M (1.4) 

the variational equation: to find v° E E p such that for all v° E H 1 ( w) N 

[ [c;^(v°)}-4(v 0 )dT = e p (v°) ( 1 . 5 ) 

J (jJ 

has a unique solution v° in E p . 

(11) Assume that uj is connected and that 70 is a subset of 7 with strictly 
positive Hn- 2 measure. Given i p E (E p 0 ) r , that is there exists c > 0 
such that for all v° E H* q (w) n 

r(,°)i<ci| £ f( V °)iu 2M ( 1 . 6 ) 

the variational equation: to find v° G E P o such that for all v° G H* g (u) N 

[ [C^e p (v°)}-e?(v°)dT = £ p (v 0 ) (1.7) 

J UJ 

has a unique solution v° in E p q . 

3 REDUCED MEMBRANE SHELL EQUATION 

The membrane shell equation can be further decomposed into a system of two 
equations. For test functions v E V, that is (vr,v°) £ H 1 { (jJ ) N x L 2 {uj) ) 

-f y° L) 2 b 

3c > 0, Vt£ G H x {uj), \£ p {v° n)| < c||t>° D 2 b\\ L ^ u ) < c' ||v°||jr, 2 ( w) 

=> 3 f p G L 2 {u>) such that t p (i>° n) = f w f p v° dT 

It will be convenient to define the function 

P( def r f p (X)/\\D 2 b(X)\\, if \\D 2 b(X)\\ t 0 

Jn [ 1 - ( 0, if \\D 2 b(X)\\ = 0 



By construction f p ||D 2 &|| E L 2 (uj). 
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Denote by H^(uj) n (resp. the subspace {u E f/’ 1 (u;) jV (resp. 

H * q (lj) n ) : v • n — 0} of tangential vectors. The decomposition yields the two 
equations 



Vfp E 



[c;. 



r,o\ 



J4(v 

Su 



-D 2 b = f p = f p \\D 2 b\\ 
C e p £ r(v 0 )] - ef ( v r) ^ — ^ P i v r) 



(1.8) 



where by condition (1.4) on i p 

3c > 0, Vt>p E H}{u) N , \i P {v^)\ < c||ef( v r)IUa (w) . 

In the case of the plate ( D 2 b — 0), e p (u°) = £p(up) + v ° D 2 b = £p(up) and 
there is only the variational equation 



V^r ^ , f w \P e p e r (^r)] " £ r ( v r) ^ ^ P ( v r ) 

which completely specifies v p E H}{uj) n / ker£p (resp. H* ot (w) N ) and v ° is 
arbitrary. There is a generalization of this result without adding new conditions 
on D 2 b. The second equation (1.8) specifies the tangential part Up of u° up to 
an element of some appropriate equivalence class providing a natural decom- 
position of the membrane shell equation into an equation for the equivalence 
class of i)p and an equation for again modulo another equivalence class. In 
the case of the plate the corresponding equivalence class for u° is so big that 
there is no information on u° and we have uniqueness for Up in the case of 
homogeneous Dirichlet boundary conditions on a part of the boundary. 

Denote by [v]e the equivalence class of v in E p (resp. E P Q ). Let 



def 



uo = {# E oj : D 2 b(x) — 0} and u>+ = uj\u}q. 
For v E V (resp. V l0 ) define the function 



( \ def 
tt 5 (u) = 



C~pD 2 b - D 2 b ’ 



m LJ. j. 

in ojq 



(1.9) 



Using the identity £p(u) = £p(ur) + v n D 2 b ) it is easy to verify that for all 
v E V (resp. V l0 ) 







C:}t?(v)~P 2 b 
C~pD 2 b - D 2 b 



D 2 b, 



in w+ 
in wq 



( 1 . 10 ) 



For each fr € Hl(u>) N (resp. H^ ot (u) N ) define the tensor 



e?(v r )^e^sM) = 



sf( v r): 



C^livryp^b 2 
C~pD 2 b - D 2 b ’ 



mW+ ( 1 . 11 ) 

in wo 



the quotient space 

V p * HIM" /**$ (resp. V p = f / kerif) 



( 1 . 12 ) 




